WEAK STABILITY OF LAGRANGIAN SOLUTIONS TO THE 
SEMIGEOSTROPHIC EQUATIONS 



JOSIANE C. O. FARIA, MILTON C. LOPES FILHO, 
AND HELENA J. NUSSENZVEIG LOPES 



Abstract. In [5], CuUen and Feldman proved existence of Lagrangian solu- 
tions for the semigeostrophic system in physical variables with initial potential 
vorticity in , p > 1. Here, we show that a subsequence of the Lagrangian 
solutions corresponding to a strongly convergent sequence of initial potential 
vorticities in converges strongly in L', q < oo, to a Lagrangian solution, 
in particular extending the existence result of CuUen and Feldman to the case 
p = 1. We also present a counterexample for Lagrangian solutions correspond- 
ing to a sequence of initial potential vorticities converging in BM. The ana- 
lytical tools used include techniques from optimal transportation, Ambrosio's 
results on transport by BV vector fields, and Orlicz spaces. 

1. Semigeostrophic equations in physical and in dual variables 

Semigeostrophic equations are simplified models for large-scale geophysical flows. 
These systems were introduced by Hoskins in [12] , as part of a family of models for 
geophysical flows under approximate geostrophic balance, i.e. where Coriolis forces 
and horizontal gradients of pressure nearly balance. We refer the reader to [7] for 
a thorough account of semigeostrophy from the physical point of view. 

In the present work, we are concerned with two versions of the semigeostrophic 
system - the incompressible 3D system and the shallow water system, both in a 
bounded domain in R'^ with constant Coriolis force. We will first focus the dis- 
cussion on the incompressible case, leaving the shallow water system to Section 4. 
The incompressible semigeostrophic equation, or SG equation, has a rich mathe- 
matical structure, closely related with optimal transport theory. Written in dual 
variables, it can be interpreted as a fully nonlinear active scalar equation, where 
the transported scalar, called potential vorticity, generates the transporting velocity 
by means of a Monge- Ampere equation. Several results are available for the semi- 
geostrophic system in dual form, beginning with the work of J.-D. Benamou and 
Y. Brenier, [5], on existence of weak solutions, and including [21 UHl [HI US] ■ In [13] j 
G. Loeper proved existence of weak solutions for potential vorticities in the space 
of Radon measures. Obtaining a solution in physical variables from weak solutions 
of the dual form is both delicate and physically relevant. For potential vorticities 
ui LP , p > 1, this problem was partially solved by M. CuUen and M. Feldman in [8], 
with the introduction of Lagrangian solutions to the system in physical variables. 

This article's main concern is the weak stability of Lagrangian solutions with 
respect to perturbations in the initial potential vorticity, complementing the work of 
CuUen and Feldman. Our main result is the weak compactness in L^, for any p < oo, 
of sequences of Lagrangian solutions, obtained from sequences of initial potential 
vorticities converging strongly in i^. In addition, we also include analysis of the 
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compressible shallow water case, as in [8]. Finally, we present a counterexample 
for the extension of our main result to initial potential vorticities in the space of 
measures. 

Let 17 C M"^ be open and bounded. The 3D incompressible semigeostrophic equa- 
tions in physical variables are a system of four equations with three components of 
velocity u = (ui, U2, U3) and the pressure p as unknown. Before we write down the 
system, we introduce the geostrophic velocity = (vf, u|, 0) as = (— 92P, 0) 
and the density p by assuming vertical hydrostatic balance p = —d^p. The semi- 
geostrophic equations have the form 

, ) + {V2, -vf ) = (m2, -Ml) 

Dtp = 
div u — 0, 

where Dt — dt + u ■ y is the material derivative. The unknowns are functions of 
(x, t) G i7 X [0, T) with initial and boundary data given by 

, , u-i^^O on dn X [0, T) 

pix,0) = po{x) in n, 

where v is the unit outer normal to dfl. Taken together, these four equations 
make up a rather odd- looking system of PDEs. We have three transport equations 
for the components of Vp, in which u enters only algebraically, together with the 
divergence-free condition. That this system turns out to be solvable only becomes 
apparent after a change of variables, which expresses its dual formulation. 

In order to present the dual variable formulation of this system, we consider the 
modified pressure P, given by 

P{x,t) =p{x,t) + ^{xl+xl) 
and we rewrite the semigeostrophic system as 

DtX = J{X - x) 



(3) 



where 



div u — 

u-v^O on dnx[0,T) 
P(x,0) = Po{x) in n, 



J= I 1 


The dual formulation is obtained by switching dependent and independent vari- 
ables in the system above. More precisely, we assume that P and u are solutions of 
the system above, and, in addition, P{-,t) is convex and smooth for all t. A con- 
sequence of the convexity of P is that VP(-, t) is a diffcomorphism between and 
some subset of M^. We introduce a = VPjJxn, where the sharp indicates measure 
pushforward and xn denotes the Lebesgue measure in fl. To be precise, if fli and 
^2 are subsets of M", fi is a measure on fli and is a measure on fl2 and X maps 
r^i to a subset of fl2 then the notation Xjl/i = means that: 

fiy)diy^ f I{x{x))dii, 
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for any / G (^^(f^a). 

The measure a is called potential vorticity. We also introduce P* the Legendre 
transform of P, i.e. 

P*{X,t) = sup{x-X-P(x,t)}. 
The potential vorticity a satisfies the following system of equations: 

{dta + y ■ (Ua) = x [0, T) 

U{X,t) = J[X-VP*{X,t)] 
VP{;t)Un = a{;t) 
a(X,0) ao(X) a.e. X G M^. 

From the definition of the pushforward measure, we can see that the statement 
VP(-,t)ttxn — a{-,t) amounts to a weak form of the equation det (D^P*) = a, 
with the condition that the image of VP* is 51. This observation shows that ([4]) is 
an active scalar transport equation, where the transporting velocity is determined 
from the transported scalar by means of a Monge- Ampere equation. The derivation 
of the dual system from the physical system is a standard calculation, and it can be 
found, for example, in The key hypothesis for the validity of this derivation is 
the convexity of P, something which is preserved under semigeostrophic evolution, 
see [3]. 

The following result concerns existence of weak solutions for the semigeostrophic 
system in dual variables. Let fl he a bounded domain in M.^, which we assume is 
contained in the ball of radius > 0, centered at the origin. Let T > 0. 

We fix the bounded smooth domain 17 in physical space, the radius 5* > and 
the time horizon T > throughout the remainder of this paper. 

We require the notation and terminology associated to Orlicz spaces in the state- 
ment below. We will briefly account for the basic theory of Orlicz spaces in Section 
2. 

Theorem 1.1. [31 [51 IE' Let Pq — Po{x) be a bounded, convex function in fl, and 
let ao -DPo#xn- Suppose that G L^(M^) for some q > 1, and it is compactly 
supported. Let Rq be such that the support of aQ is contained in the ball B{0,Rq) 
and set R{T) — Roe"^ + {e^ — 1)S . Then, for any t > 0, there exist functions 
a = a{X,t) G L°°([0,T),L«(R3))^ p ^ P{X,t) G L°°{[0,T),W^'°°{n)), and an 
N -function A such that 

(i) 

(5) supp{a{-,t)) C P(0,i?(r)); Vt G [0,T), 

(ii) For each < t < T , P(-, t) is convex, and a, P satisfy 

aeC{[Q,T),LA{B{0,R{T)))) and P e C{[0,T),W^^^{n)), 

for any real number r G [l,oo), where La is the Orlicz space associated 
with the N -function A. 

(iii) P*{-,t) is convex, pointwise in time, locally bounded in space-time and 

VP* G L°°([0, T),Ea' (P(0, R{T)))) n C([0, T), L'-(P(0, P(T)))) 
for any r G [1, oo), where Ea' is the dual of La- Moreover, 

(6) llvp*(-,t)|U.o(R3) < 5 v<G[o,r), 



4 



HELENA J. NUSSENZVEIG LOPES ET AL. 



(iv) (a, P, P*) satisfy where the evolution equation and the initial data for 
a are understood in the weak sense, i.e., for each (j> G C^{M.'^ x [0,T)) 



[dt(j){X, t) + U{X, t) ■ V0(X, t)]a{X, t)dXdt+ 



(7) 



R3x[0,T) 



+ / ao{X)(f>{X,0)dX = 0. 



One would like to find solutions for the SG system in physical space, which 
means, solutions to ([3]). Our point of departure is the work of CuUen and Feldman 
in [8]. In that article, CuUen and Feldman pointed out that: 

• concerning the Eulerian form of system ([3]) - a distributional formulation 
of this system requires making sense of products of components of u with 
first derivatives of P; 

• the formal expression for u is given by 



• the known regularity for solutions of the dual problem has P* Lipschitz 
continuous. 

Clearly, making sense of the physical velocity u is complicated, given that D^P* 
is a measure while VP is only bounded, not to mention making sense of the product 
u ■ VP. As a consequence, seeking Eulerian solutions in physical variables is a 
difficult problem. This was the motivation for the introduction of the notion of 
Lagrangian solutions. 

Let Po e be a convex function and r e [l,oo). Let P : Q x [0, T) R 

satisfy 



Definition 1.1. (Lagrangian Solutions) The pair (P, P) is called a Lagrangian 
solution of in fl X [0, T) if 

(i) P(x, 0) = X, P{x, 0) = Po(x) for almost all x € fl, 

(ii) for all < t < T the mapping Ft = F(-,t) . CI —> preserves Lebesgue 
measure, i.e. Pf#xo = Xn, 

(iii) There exists a Borelian map F* : Q x [0,T) ^ fl .such that, for all t G 
(0,r), the map F^ = F*{-,t) = f2 — > 17 preserves Lebesgue measure, (i.e. 
-Pt*#Xf2 = Xn), and satisfies Ft o Ft{x) = x and Ft o Ft{x) = x for almost 
all X £0,, 

(iv) The function 



u{x,t) = dtWP*{\/P{x,t),t) + D^P*{VP{x,t),t)[J{VP{x,t) -x)]; 




(11) 




(12) 
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in the following sense: for any ip ^ C^{fl x [0, T)), 

[Z{x, t) ■ dtifix, t) + J{Z{x, t) - F{x, t)) ■ (p{x, t)]dxdt+ 



(]^3) Jnx[o,T) 

\7Pq{x) ■ (p{x,0)dx = 0. 



Given a Lagrangian solution {P,F), the map F{-,t) is called a Lagrangian flow 
in physical space, for each t e [0, T). 

Next we give the precise statement of existence of Lagrangian solutions. 

Theorem 1.2. Let Pq be convex and bounded in B(0, S). Assume that Pq satisfies 

(14) DPo#xn e L''{m.^) 

for some q > 1, and that DPoffxii compactly supported. Then there exists a 
Lagrangian solution {P, F) of (0) in Q x [0,T), for which ^-(W^ are satisfied 
for any r G [l,oo). Moreover, the function Z = Z{x,t), defined by satisfies 
Z{x, ■) G M^^'°°([0, T)) for almost all x G and iPJ\) is also satisfied in the following 
sense: 

, . dtZ{x,t) ^ J{Z{x,t) ~ F{x,t)), for {x,t) enx {0,T), -a.e. 
^ ' Zix, 0) = VPo(x), for xen, - a.e. 

The case g > 1 is the main result in |B] . It was observed by one of the authors of 
the present paper that the result in [8] can be extended to g = 1 using the technique 
from [M] . The proof of the case g = 1 of Theorem 11.21 is embedded in the current 
work. 

Let us briefly examine the construction underlying the proof of Theorem 11.21 
Under the hypotheses of Theorem 1.2, from the solution of the dual problem, and 
using the Polar Factorization Theorem (see [3]), one obtains P. Given P, it can 
be shown that the following expression gives rise to a Lagrangian flow in physical 
space: 

(16) F = F{x, t) = VP* o $t o VPo{x), 

where, for each t, is the Lagrangian flow in dual space, obtained using Ambrosio's 
theorem as follows. 

Consider the transport equation 

dta + U.Va = 0, 

which is equivalent to the first equation in ([4]), since div [/ = 0. From the regularity 
of P* we have 

Ui^LToM^ X [o,T)), c/eL°°([o,r),PKoc) 

One uses Ambrosio's theorem to obtain the Lagrangian fiow associated with the 
transport equation above. To do so one must modify the velocity U near infinity 
without affecting the solution a; this can be achieved since a has compact support 
in X [0,r). We have 



supp a C B{0,R{T)) x [0,T], 

where R{T) was introduced in Theorem ll.il (i). 

We introduce a modified velocity U as follows: choose g G C°°(R) such that 

(17) = 1 in {|s| < i?(r)}, e = in {|sl > P(T) + 1}, < g < 1 in R, 
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and define, for X e M'', 

(18) H{X) = (e(|Xi|)Xi, q{\X2\)X2, q{\X^\)X^). 
The modified velocity U is then given by 

(19) U{X,t)^J[H{X)-VP*{X,t)l 
and, therefore U satisfies 

U e i°°(]R3 X [0,T)) 

(20) e L°°([o,r),Bv;„J 

div tj{-,t) = in R^, for all t e [0,T). 

Furthermore, 

(21) U^U in B{0,R{T)), 

(22) ||[/(-,i)|Uoo(R3) < 5 + i?(r) + l forall [0,r). 
Therefore, using the results of L. Ambrosio in [2\, we have 

Proposition 1.1. Let U be given by 119\) . There exists a unique locally bounded 
and Borel measurable mapping $ : x [0, T) — > satisfying: 

(i) $(X, •) e H/i'°°([0,T)) /or aZmosf X £ M^. 

(ii) $(X, 0) = X, X G £3 _ a.e.; 

(iii) For almost all {X,t) G ^ (o,r) 

(23) dt'^{X,t) = U{<^{X,t),t)- 

(iv) <i>(-,t) : M'' ^ M'^ preserves Lebesgue measure in for all t G [0,T). 

(v) We have 

(24) ^iX,t) C B{0,R{T)) for almost all {X,t) G VPo(f^) x [0,T). 
/n particular, 

(25) 9t$(X, t) = C/(«>(X, t), i) /or almost all {X, t) G VPo(f^) x [0, T). 

(vi) There exists a Borel mapping $* : M'^ x [0, T) — > 'M? such that, for all 
t G (0,T) f/ie map : M'^ — > preserves Lebesgue measure in M"^, and 
satisfies '^l o $4(2:) = a; and $t o <l>j (x) = x, for almost all x , 

(vii) Under the conditions of Theorem 1.2, if {ct,P) is a weak solution of 
and if (i)- (vi) hold, then for any t G [0,T], 

(26) at = $t#ao. 
Moreover, for any t G [0,T], 

(27) at{x) = ao(^>t (a;)) for almost all a; G M^. 

2. Orlicz spaces 

In what follows, we will collect definitions and a few results concerning Orlicz 
spaces. For details, see [T] and [13]. 

Consider a : [0, 00) — + [0, 00) with the following properties: 

(i) a(0) = 0, a{t) > if i > and lim a{t) 00, 

t — 'OC 

(ii) a is non-decreasing, 

(iii) a is right-continuous. 
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The function A, defined on [0, oo) by taking 



A{t) = / a(r)dT, 







is called an iV-function. 

We note that iV- functions are continuous on [0, cxd), convex and strictly increas- 
ing. 

An A^-function is said to be A-regular if there exists a positive constant C and 
to>0 such that A{2t) < CA{t), \/t>to. 

Let be a domain in i?", and A an A^- function. The Orlicz space La{^) is 
the linear closure of the set of functions u G Lj^^{n) such that ^(|w|) is integrable. 
These are Banach spaces with norm given by 

|u|U = inf <i A: > 0, f A ( ^-^^] dx < 1 



Jn \ 1^ / 

They generalize the Lebesgue spaces Lp{^1), which are Orlicz spaces, with A^- 
function A{t) = tP. We denote by Eji{il) the closure, with respect to || • ||^, of 
the set of smooth, compactly supported functions in Q. For every A'^-function A, 
we have that Ea{^) is separable. In general, La{^) and Ea{^) are distinct, and 
L^(f2) is not separable. However, when A is A-regular, La = Ea- 

Let A be an A^-function. Its Legendre transform A* is given by 

A* = A* (s) = maxlst - A(t)}. 

t>o 

One can verify that A* is also an A^-function and that A** — A. 

Finally, the following classical results will be relevant in our analysis. 

Theorem 2.1. [IJ The dual of Ea{^) is LA*{n). 

Lemma 2.1. [6^ Let il be a bounded domain in and f, G L^iVl), for all k. 
If f'^ ^ f strongly in L^ then there exists a A-regular N -function A such that {/'"'} 
and f are uniformly bounded in i^(f2). 

Lemma 2.2. |14| Let {u„} be a sequence of functions uniformly bounded in L°°(fl). 
If Un —>■ u in L^(n), then Un ^ u in La{^), for any N -function A. 

3. Weak stability of the semigeostrophic equations in physical space 

We now turn to the main objective of this article. We consider a sequence of 
initial potential vorticities, converging strongly in L^ to a given limit vorticity. 
We would like to understand the convergence properties of the corresponding La- 
grangian solutions in physical space. Our motivation for considering this problem 
was, originally, to try to extend Cullen and Feldman's construction to solutions of 
the semigeostrophic equations with measures as potential vorticities. To do so we 
intended to approximate such solutions by smoother ones and, hence, we needed 
to understand how the corresponding Lagrangian solutions behaved. As it turns 
out this approach to construct Lagrangian solutions for measure- valued potential 
vorticities does not work; this will be made clear by means of a counterexample, in 
Section 5. Instead, we have established a weak stability, or continuity property, of 
Lagrangian solutions with respect to integrable perturbations of an L^ initial po- 
tential vorticity. Weak stability of weak solutions of the semigeostrophic equations, 
in dual formulation, has already been established by G. Loeper in [13] . 
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Throughout the remainder of this section we fix ao G i^(M'^) together with a 
sequence {oq} C i^(M'^) such that ao strongly in L^(M'^). In addition, we 

assume that ao and the sequence {ckg} are compactly supported, with supports 
contained in a ball B{0,Ro). 

Using Lemma [2.11 as in [M], we have that there exists a A-regular iV- function 
A such that ag, ao are uniformly bounded in 

Let a" = a"(?/,t) be a weak solution of the semigeostrophic equations in dual 
formulation with initial potential vorticity ag . Consider the corresponding modified 
pressures P", Pq, defined in the physical space fl. Denote by <I>" the Lagrangian 
flow in dual space given in Proposition 11.11 Finally, consider the corresponding 
Lagrangian flows in physical space, 

(28) F," := ViPn* o $r o VPg", 

as obtained in Theorem II. 2[ equation 

In the proof of Theorem ll.li case q = 1, (see [l^), it was shown that a sub- 
sequence of {a"}, {P"}, {(P")*} exists, which we do not re-label, together with 
a weak solution a = a{x,t) of the semigeostrophic equations in dual formulation 
(with initial potential vorticity ao), such that the following hold, for each t G [0, T): 

a"(-,t) -a(-,i) u.*-L^(P(0,P(r))) 

p"(-,<) — >P{-,t) w'^^^'in) 

(29) Po"— ^^^0 w^^^n) 
(P")*(.,t)-^P*(.,<) <'i(K3) 
V(P")*(.,i) ^ VP*(-,t) {Ea')ioc{M.^)- 

In the proof of the convergence of V(P")*(-,t) to VP*(-,t) one uses Theorem 12.11 
and Lemma [121 

We fix, throughout the remainder of this section, such a subsequence. 

Above, P(-, i), Po are the modified pressures corresponding to a(-, t) and ao, and 
R{T) is given in Theorem ll.li item (i). Let $ be the Lagrangian flow in dual space 
as in Proposition 11.11 and let 

(30) Ft VP;o$toVPo, 

be the Lagrangian flow in physical space obtained in Theorem 1 1.21 see 
Our main result is the following. 

Theorem 3.1. There exists a (further) subsequence {P"''} C {P"} such that, for 
almost every < t < T , we have: 

P"'' — > Ft, strongly in L''(f2), as k oo, 

for all r G [1, oo). 

Before we present the proof of the theorem, we require the following auxiliary 
result. 

Lemma 3.1. For each R > 0, we have 

lim / sup |$"(X,t) - $(X,t)|dX = 0. 

Jb(o.r.) telo.T] 
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Proof. For each n, recall that <&" is the Lagrangian flow in dual space associ- 
ated with the vector field C/"(X,t) = J[H{X) - V(P;')*)]. Since we have that 
(P")*(-,t) ^ P*(-,i) in W^itc (K^), it follows that 

(31) U-{X,t)^U{X,t)emLlM^). 

However, such a condition is not enough to obtain the convergence of due 
to the fact that we cannot control VC/". This is required in Ambrosio's stability 
result, namely Theorem 6.5 of [5]. 

We consider, instead, a family of aproximations of C/" given by 

(32) i7"^"(X,t) - J[HiX) - {\7{Pir*vmX,t)], 

where 77™ is a standard moUifier. 
Now we have: 

jjn,r,i g C{[0,T] X R^,R^), 

SUPm l|C^"'"'llL=°(R3x[0,r],R3) < OOi 

(33) div C/"'™ = 0, 
||VC/"'™||ioc([o,T]xB(o,fl),R3) < C{n,m,R) < 00, 

ljn,rn , (jn il^^(]R3 ^ [0,T)). 

Let t) be the Lagrangian flow associated to C/"'™. It follows from The- 

orem 6.5 in [2] that 

(34) lim / sup |$"(X,i)-$"'™(X,i)M^ = 0, Vi?>0. 

™^°°JB(0,i?) te[o,T] 

Note that U'^-"'{X,t) U{X,t) in L,^o^(M^ x [0,r)) when m,n ^ 00. To see 
this it is enough to observe that, for any i? > 0, we have 

l|v(p,")* - vp;iUi(B(o,i?)) < lb7'"llLHs(o.fl))l|v(pr)* - vp;iUi(B(o,i?)) 
+11 vp; * 77" - vp;iUi(B(o,/?)) 0. 

Hence, we also have 

(35) lim / sup |$(X,t) - $"'™(X,t)|dX = 0, VP > 0. 

Js(o,fl) te[o,T] 

Given ([M)) . it is possible to choose a subsequence m — m{n) > n such that 

(36) lim / sup 1 4>"(X,t) I dX = 0, VP >0. 

"^°°JB(0,fl) te[o,T] 

We conclude, from and that 

lim / sup |$"(X,i) - $(X,i)|rfX = 0, VP > 0, 

tG[0,T] 

which concludes the proof. 

□ 

Remark 3.1. Once we take into account the expression U6]) . Proposition 1.1 (v), 
and (23), we see that we may assume in what follows that the flow ^{X,t) is 
associated with the vector field U{X,t) — J[X — \7P*{X,t)]. 
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With this lemma we are now ready to give the proof of Theorem 13. II 

Proof of Theorem \3.1\ Let us first prove our result for the case r = 1. 
We note that, for each < t < T, we have 



(37) 



Fl'{x)-Ft{x)\dx = 
Jn 

< / |V(i^")*o$'/oVPo"(a;)-V(Pt")*°*t°VPo"(a;)|rfa; 
Jn 

+ / |V(Pt")* o $i o VPo"(x) - VP; o $t o VP^{x)\dx 
Jn 

+ / |VP; o $t o VP^ (x) - VP; o $f o VPo(a;)| 
Jn 



In 



We will show that each of these integrals vanish as n tends to infinity, passing 
to subsequences as needed. 

Let us begin by considering Ii. Using that VP(5*jJxr2 — C(q we have: 

/ \v{pn* ° <fr ° vPo"(x) - v(pr)* o o vPo"(x)i dx = 

Jn 

|V(Pr)* o cO«(y) - V(i^")* o <i>t{y)\a^{y)dy 



(38) </ |V(Pr)*o<DJ'(y)-VP;o$«(y)|aJf(y)dy+ 



|VP; o - V(Pr)* o $,(y)| a^(y)d2/ 

From Proposition l.f (vii) we sec that a" = $"#aQ and, since La is a rearrange- 
ment invariant space, it follows that l|ao||L^ — for each n and for each 
t e [0,r). Therefore, 

II - |V(Pr)* o ^\y) - VP; o ^\y)\ al(y)dy 



(39) = \y{prn^) - yp:iz)\a^{z)dz 

<||v(pr)*-vp;ib,.||ariu, 

= l|v(pn*-vp;ib,.||a[jlU,"=5Po, 

where we have used ((29|) . Theorem 12.11 and the boundedness of ag in La- 
As for we have: 

I! |VP;o<J>J^(y)-V(P,To<J>,(y)|a^(y)dy 

Jf 

/ |VP;o<i>J'(y)-VP;o$,(y)|aS(y)dy+ 

+ / |VP; o $,(y) - V(Pr)* o $,(y)| a^(y)dy 

_ r2,l I t2.2 
— ^1 I -"l ■ 

Consider the integral I^'^. Since — ^ ao strongly in L^, it is easy to see that 

the proof that /^'^ tends to zero as n — * cxo reduces, by Lebesgue's dominated 
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convergence theorem, to showing that, for all t, 

(41) VP;o$J'(y)- VP;o$t(y) — >0, n->oo a.e. y G 

At this point we must pass to a further subsequence. We have, by Lemma |3.1[ that 
5":=.9"(y)= sup ^0 strongly in X;^^. 

0<t<T 

From this it follows that there exists a subsequence, {5"''}, which converges, a.e. 
y e M^, to as fc — > cxD. Hence, for every < t < T, we have 

(42) ^^'(y) - '^t{y) a.e. ?/ e as fc 00. 

Now, since P^* is convex, it follows that VP^* is almost everywhere differentiable 
(see, for instance, [11]), and hence, continuous except for a set of Lebesgue measure 
zero, say N C R^. Given (|42|l it is enough to show, therefore, that, for almost all 
2/ G R"^, VPj* is continuous at $t(?/). To see this we note that 

\{y G R^ My) e N}\ - \{y e R=^; y G $:(A^)}| = |A^| = 0, 

in view of the fact that $j preserves Lebesgue measure. Therefore, $t(2/) is a 
continuity point for VP^*, for almost all y, as desired. 
Hence, 

(43) / |VP;o$;''=(y)~VP;o$,(y)|a^'=(y)dy^O. 

The analysis of I^'^ is similar. We have that VP;(R3), V(P4")*(R3) c P(0,5), 
Vt, 71, V(Pt")* VP; strongly in L^^, hence V(P")* -> VP* strongly in P,^^^(R^ x 
R+). Thus, we may pass to a subsequence, chosen independently of t, and which 
we do not re-label, so that 

V(P"'^')* ^ VP; a.e. yGR^ 
as fc ^ 00, for almost every < i < T. 

Using this, together with the fact that <i>j is measure preserving, we may conclude 
as before that 

ViPtT" ° My) - VP; o My) ^ O, a.e. y g R^ 
a.e. < t < P. This, together with the strong convergence in of ao, and 

Lebesgue's dominated convergence theorem, yield, 

(44) / \v{p:r''oMy)-^P:°My)\ao''iy)dy^O- 

From (|43p and we have that /J — > 0, which concludes the analysis of Ji. 
Next we consider l2- Using the fact that VP^ij^Xit — '^oi have that, 

( |V(P;')* o $t o VPo"(a;) - VP; o $t o VPo"(x)| dx = 

(45) -^f^ . 

= / |V(P,")*o$,(y)-VP;o$,(y)|ao"(y)dy 

which is the same as /^'^. Hence, from ((44)) . it follows that, passing to the appro- 
priate subsequence, /2 — * 0. 

Finally, we consider the last integral, 

h = I |vp; o o vPo"(x) - vp; o $t o vPo(x)| dx. 

in 
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Now, since VP* is bounded, it is enough, by the Lebesgue dominated convergence 
theorem, to prove that 

(46) VP; o $t o VFq"" {x) - VP; o $t o VPo{x) O a.e. x e fi. 

To this end we, once more, pass to a subsequence for which VPq'' VPq a-e. 
X e il, and we do not further re-label. 

Recall that the support of ap was assumed to be contained in the ball B{0, Rq), 
for all n, and that VPo'(^) is precisely the support of Qq, hence contained in 
B{0,Ro). 

Next, we note that, passing to subsequences as needed, $f o VPq"*" ^ o VPq 
a.e. X G ri. We will show this by proving the convergence of <i>f o VPq"* ^$40 VPq 
in and passing to a subsequence which converges a.e. x fl. 

By Lusin's theorem we have that $t coincides with a continuous function up to 
a set of arbitrarily small Lebesgue measure. More precisely, let £ > and consider 

e C°{B{0,Ro)) andP= C P(0,Po) such that $t = f outside of P= and |P=| < e. 
Since <i>t is bounded, for each t, we may assume that is also bounded, uniformly 
in e. We use the fact that VP|^''"#Xf2 = chq'' , and the analogous fact for Pq, to 
estimate: 

limsup [ |$t o VPq"" - $t o VPoI dx < limsup / |($t - /^) o VP^" \ dx 
+ limsup / \f o VPo'"^' - o VPol dx + /" |(/^ - $() o VPol dx 

rifc^cxD Jo Jo 

= limsup/" \^t~ r \ da'^i" +limsvLp [ \roVP^''^f oVPo\dx+ [ \f-^t\dao 

rifc— ♦00 JiJe rifc— ►00 Jo J_B= 

<2||$,-/-||i^(|a^'^(i?-)| + |ao(i?")|) + limsup / 1/^ o VPq"'' - ^ o VPo| dx. 

rifc^oo Jo 

The first term can be made arbitrarily small since {a^} is uniformly integrable, 
while the second term vanishes because is continuous. 

We have shown that <I>t o S/P^'' — > <i)j o VPq a.e. x S fi, passing to a further 
subsequence if needed. Next, recall that VP^* is continuous in M.^ \ N, so that, 
passing to the subsequence above, to obtain it is enough to show that 

\{xen; $toVPo(x) eiV}| =0. 

Recall that preserves Lebesgue measure, so that |$( (iV)| = \N\ — 0. With this 
we obtain, using again that VPo#xo = ceo: 

\{x en;^to VPo(x) e iV}| = \{x e n; VPo(x) e $*(iV)}| 
(^^y^ = ^ X*:(iV) ° VPo(x) dx 

= / ao{y)dy = 0, 

and therefore, I3 0. 

This establishes our result if ?■ = 1. 

Now, given that V(P")* is uniformly bounded, we obtain the convergence in 
L'"(ri), 1 < r < cx), by interpolation. This concludes the proof. 
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4. Weak stability for the shallow water case 

The shallow water version of the semigeostrophic equations can be written as an 
equation for h — h{x,t), x = {xi,X2) G C M^, t e [0,r), and v — (wi,W2)- Here, 
h is the height of fluid above H. and v is the velocity. We denote Dt ~ dt + v ■ V 
and we set 

P = Pix,t) = h{x,t) + ^\x\^ and VP. 
With this notation, the shallow water SG equations take the form: 



(48) 



DtX = {X -x)^, nx{0,T), 

dth+ div{hv) = 0, nx (0,T), 

viy = 0, on dflx [0,T), 

h{x, 0) = ho{x), in fl. 

Here, {a,b)-^ = (—6, a). 

In dual variables, this problem can be written as 



dta + V ■ (Ua) = 0, x (0, T) 

VPt^ht=at; ie(0,T) 

(49) ) U{X,t) = [X-VP*iX,t)]^ M2x[0,T) 
P*{X, t) sup{x ■ X - P{x,t)}, M2x[0,T) 

a(-,0) ^ao = VPo#ho. 

For the modeling background concerning this system, see [H [9] . A weak solution 
for this system was obtained by M. CuUen and W. Gangbo, see [9], in the casep > 1, 
and their existence result is similar to Theorem 1.1. In [5], Cullen and Feldman 
also proved existence of Lagrangian solutions in physical space for the system p8)) 
for p > I. The existence results, both for weak solutions in dual variables (from 
[5]) and for Lagrangian solutions in physical variables (from 8J), can be extended 
to p = 1. The proof for weak solutions is an easy adaptation of the work in T^, 
whereas the proof for Lagrangian solutions is, as before, embedded in what follows. 

We are interested in Lagrangian solutions {P, F), where _F : x [0, T) ^ is a 
Lagrangian flow associated with v, and P is obtained from a weak solution in dual 
variables. However, for the shallow water case, the vector field v is not divergence- 
free. Nevertheless, the transport equation dth + div{hv) — holds. Therefore, if 
_F is a Lagrangian flow associated with v, the solutions h of this equation satisfy 
Ft4t^ho — ht, Vt S [0,r). This property replaces the fact that F preserves the 
Lebesgue measure in the incompressible case. 

Let C be open and bounded and let T > 0. Let Pg = -Po(^) be a convex, 

bounded function in such that ^.0(2;) = Pq{x) — 2 1^1^ > in f2. Let r e [l,oo) 
and F : 17 X [0, T) ^ M be such that 

(50) p e i°°([o, T), w^'°^{n)) n c([o, r), w^^^'in)) 

(51) -P(-,i) is convex in Q for each te[0,T). 

Let h{x,t) — P{x,t) — Let F : x [0, T) be a Borel map satisfying 



(52) 



F e C{[0,T),L'-{n,hodx)). 
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Definition 4.1. (Lagrangian Solutions) The pair (P, F) is called a weak Lagrangian 
solution of ^ in VI X [0, T) if 

(i) F{x,Q) = X, ho-a.e. inQ, P(x,0) = Po{x) a.e. inVt, 

(ii) for every t > Q the map Ft = F{-, t) : Q ^ n is such that Ft^ho ~ ht, 

(iii) There exists a Borel map F* : Q. x [Q,T) ^ such that, for each t e (0, T) 
we have F^ = F*{-,t) = Q satisfies F^^ht — ho, and Ft o F^{x) — 
X ht — a.e. in fl and F^ o Ft{x) = x h^ — a.e. in il, 

(iv) The function 

(53) Z{x,t)^VP{Ft{x),t) 

is a weak solution of 

, - dtZ{x^t) = [Z{x,t) — F{x^t)]^ ^ on supp h^ in J7 x [0, T) 
Z(a;, 0) = VPo(2;), on supp h^ in ft, 

in the following sense: for every ip G C^(il x [0, T)), 

[Z{x,t) ■ dtip{x,t) + {Z{x,t) - F{x,t))^ ■ ip{x,t)]hQ{x)dxdt+ 



(55) Jn^T) 

VPoix) ■(p{x,0)hoix)dx = 0. 
The map Ft is a Lagrangian flow in physical space. 

With this definition in place we give the precise statement of the existence of 
Lagrangian solutions in the shallow water case. 

Theorem 4.1. Let C be open and bounded, and assume that U, G B , where 
B is the open ball B{0, S). Let hQ{x) > be such that Pq = Po{x) = ho{x) + ^\x\'^ 
is a convex, bounded function in B , and assume that 

(56) DPo^ho e L«(VPo(^^)) 

for some q > 1. Then, for each T > 0, there exists a Lagrangian solution {P, F) 
of |^($[ ) m il X [0,r), where (50)~(52) hold for all r G [l,oo). Furthermore, the 
function Z — Z{x,t) defined in (EH) satisfies Z{x, ) G W^'^(\Q,T)) hg-almost 
everywhere in fl, and ^54^ is also satisfied in the following sense 

. . dtZ{x,t) = {Z{x,t) - F{x,t))^, haC^ X a.e. in n x {0,T), 

^ ' Z{x,0) ^yPoix), hoC^ a.e. in n. 

As before, one obtains P from the dual problem; CuUen and Feldman showed 
that F(-,t), given by the expression below, is a Lagrangian flow in physical space: 

Fix,t) ^VPt* o^toVPoix), 

where, for each t, <I>t(X) is the Lagrangian flow in dual space associated to the 
vector field U{X,t) — [H{X) — VPtiX)]^, whose construction in arises in the 
same manner as for the incompressible case. 

Let us now address the stability of Lagrangian solutions. Consider ao,aQ G 
L^(R^) with ag — > ag in L^ and with supports contained in a single ball B{0,Ro). 
Let a" = a"'{x,t) be weak solutions in dual variables with initial data aJJ. As 
before, there exists an Orlicz space La, with A-regular TV-function A, such that 
ao is uniformly bounded in La{R^). 

Let a" = a"{x,t) be a weak solution of (|49p with initial data aJJ and let /i" be 
the corresponding height and P" be the corresponding modified pressure. It can 
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be easily deduced, from the proofs of Lemma 3.6, Lemma 4.3 and Theorem 4.4 of 
[9] that, since ctg — > ao strongly in L^, there exists a subsequence such that: 

^ wc&k- *BM, 

h'^i;t)—.h{;t) in L-(17), 

^ ' (Pp)* > P; strongly in W^^l , 

V{Pp)* VP; weak - 

for each < i < T, with a, ft,, P* a weak solution of the semigeostrophic shallow 
water equations. 

We use to denote the Lagrangian flow in the dual space associated to C/", 
and we denote by i^" := V(P")* o $" o VP^ the corresponding Lagrangian flow 
in physical space. Accordingly, let $ denote the Lagrangian flow in dual variables, 
associated to the limit velocity U and let Ft := VP^* o <i)j o VPq be the corresponding 
Lagrangian flow in physical space. 

We note that the result in Lemma 3.1 remains valid in the present case. 

Theorem 4.2. There exists a subsequence {P"''} C {P"} such that, for almost 
every t G [0,r), we have 

(59) lim / IF^" (x) - Ft{x)\'' ho{x)dx = 0; 

for any r G [1, oo). 

Proof. Since F" is bounded uniformly in L°°([0,r) x fi), since Hq Hq uniformly, 
and since f2 is bounded, it is clearly enough to prove that: 



(60) lim / \Fl\x)~Ft{x)\''h'f^{x)dx = 0. 



(61) 



To show ()60p we note that, as in the incompressible case, we need only analyze 
the case r = 1, as r > 1 follows by interpolation. We have: 

\F^ix)-Ft{x)\h^,{x)dx < 

< 1^ \w{pn* o <i>- o vp-(x) - v{p-r o o vp-{x)\h^{x)dx 

+ I |V(P,")* o $t o VPo"(x) - V(P;) o $t o VP,'lix)\h^{x)dx 
+ J |VP; o $t o VPo"(^) - VP; o $t o VPo(x)|/i^(a;)dx| 

= Ui +/2 +/3}- 

The analysis of each of these integrals follows closely the analysis performed on 
the analogous integrals in Proposition 3.1, once we use the facts that VPq^Hq = 

5. An EXAMPLE IN THE SPACE OF MEASURES 

The purpose of this section is to describe a counterexample for Theorem 13.11 
for potential vorticities which are not absolutely continuous with respect to the 
Lebesgue measure. 

For the discussion in this section we will ignore the vertical variable in the in- 
compressible SG equations; the argument we will present can be easily adapted to 
accomodate the third direction. 
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We fix the physical space to be the planar disk = B{Q, 1). Let zq = (f , 0) and 

set 

ao = 7r(5z(,, 

where (5p denotes the Dirac measure at P. Let Q!(t) = 7r(5j,(t), with = (cos<,sint). 
ft can be checked that a is a weak solution of (jH), in the sense of [13]. Next, observe 
that the unique (up to a constant) convex potential for the optimal transport map 
between xo. and a{t) is given hy P = P{x,t) = z{t) ■ x. Its Legendre transform 
is P*iy,t) = \\y- z{t)\\ and, consequently, VP*{y,t) = nfff^y^- The Lagrangian 
flow in dual space, restricted to the support of ap, is precisely zq i-^ z{t). The 
Lagrangian flow in physical space cannot be computed by since $i o VPo(') 
is identically equal to z{t), where VP* is not defined. 

One can use approximations as a strategy to circumvent the difficulty described 
above; as we will show, this does not work. 

Proposition 5.1. Let z{t) — (cost, sin t) and set 

a = ■^XB(z(t),e)- 

Then is an exact weak solution of the semigeostrophic equations in dual variables 
@ with initial potential vorticity a^(-,0). 

Proof. Let us first establish the relation between a potential vorticity of the form 
and the corresponding velocity [/^. To this end, we fix ^ S and we consider 



1 



a 



The optimal transport map between xn and is given by VP , where the convex 
potential P is, up to a constant, 

P' = p'{x) = z-x + e'^-^, 

and hence VP'^{x) = z + ex. Indeed, it can be easily verified that VP^#xo = 
and P is convex, so that the uniqueness part of Brenier's Polar Factorization 
Theorem, see [1], may be applied. The Legendre transform of P is 

\\y-n' 



2^ , ifyeB{z,e) 
ify^B{z,e) 



Therefore, we find that 



(62) V(P^)*=V(P')*(t/) 



y- z 



y- z 



if y G B{z,s) 



iiy^B{z,e). 



\\y - z 

For each fixed t, we have that is of the form with z = z{t). Therefore the 
corresponding semigeostrophic velocity U^, in dual variables, is given by 

U'^U'{y,t) = (y-W{Pn*{y))^, 
where \7{P^)* is given by the expression in ([62]) with z = z{t). 
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Consider yo G B{zo,e). Let y = y{t) be the solution of 

y'^U%y,t), 
2/(0) = ya- 



As long as y{t) 6 B{z{t), e) we see that 

2/ = — 



z(t) 

e 



We also have 

z' = z^. 

Thus, subtracting these two equations, we deduce that 

{y - z)' = - — -{y - z)^, 
e 

[ iy~zm = yo-zoeBiO,e). 

Therefore y — z rotates around the origin at the rate (e — l)/e hence, in particular, 
y{t) rotates around z{t) and never leaves B{z{t),e). We have shown that the flow 
of maps B{zQ,e) to B{z{t),e) through a rigid rotation. 

This implies that a' is a weak solution of the transport equation dta'^ + U' -Va^ = 
0, as desired. 



□ 



Remark 5.1. Note that Q!^(-,<) ^ a(t) weak-* BA4, in accordance with '13]. 



Remark 5.2. From the proof above we obtain an explicit expression for the La- 
grangian flow in dual variables for Lagrangian markers inside B(zQ,e), namely: 

cos(^t) -sm{^t) 



(63) 



$^ = ^Uyo) - z{t) + 



sm 



(yo - zo). 



-t) cos(i^t) 

Next we compute the Lagrangian flow in physical space associated to a"^ using 
expression (fTH)) . Let x £ fl and note that VPo(x) — zq + ex £ B{zo,e). Hence 
we may use the Lagrangian map (|63p . together with the expression in (j62[) with 
'z — z(t), to obtain: 

(^i) -sin(^t) 



(64) 



F,'^F,^x) = 



cos I 



.n(^t) cos(^t) 

In other words, as e — > 0, describes a rotation around the origin in physical 
space with arbitrarily large angular velocity. In short, a concentrated vortex in dual 
space corresponds to a Lagrangian fast eddy in physical space, but concentrating 
the dual space vortex into a point produces an unphysical eddy which rotates 
at infinite speed. This shows that it is impossible to extend the weak stability 
theory we developed here in to the full space of measures, while keeping the 
strong convergence of sequences of Lagrangian flows as a conclusion. There are two 
possibilities for further work in this direction. One is to develop a theory of weak 
convergence of Lagrangian flows associated with converging sequences of potential 
vorticities in the space of measures and another is to try to extend the theory to 
spaces of continuous measures, considering that Diracs in dual space are associated 
with unphysical infinite velocity eddies and are, therefore, unphysical themselves. 
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but perhaps other measures, such as potential vortex sheets, may be associated 
with meaningful flows. 

We conclude with the following remark. We established the convergence of La- 
grangian flows a.e. in time, in space. However, this may not be optimal, and 
this leads to an interesting line of investigation. It was pointed out, by Brenier and 
Gangbo in [5 , that the topology induced by L''-convergence in the space of dif- 
feomorphisms is not very satisfactory. One may investigate, for instance, whether 
the convergence of Lagrangian flows can be improved for potential vorticities in 
Holder spaces, using the regularity theory for optimal transport developed by Ma, 
Trudinger and Wang in [15] . 
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